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We consider the problem of bound states in strongly anisotropic ferromagnetic impurities in a 
superconductor, motivated by recent experiments that claim to observe Majorana modes at the 
ends of ferromagnetic wires on a superconducting substrate [S. Nadj-Perge et al, Science 346 , 602 
(2014)]. Generalizing the successful theory of bound states of spherically symmetric impurities, we 
consider a wire-like potential using both analytical and numerical approaches. We find that away 
from the ends of the wire the bound states form bands with pronounced van Hove singularities, 
giving rise to subgap peaks in the local density of states. For sufficiently strong magnetization of 
the wire, we show that this process generically produces a sharp peak at zero energy in the local 
density of states near the ends of the wire. This zero-energy peak has qualitative similarities to the 
claimed signature of a Majorana mode observed in the aforementioned experiment. 

PACS numbers: 74.25.Ha,73.20.Hb 


Introduction. The antagonistic relationship between 
superconductivity and ferromagnetism manifests itself 
in many diverse ways. For example, it is well known 
that a magnetic impurity in a superconductor supports 
bound states with energy within the superconducting 
gap [1-3]. These Yu-Shiba-Rusinov (YSR) states can 
be directly imaged using scanning tunneling microscopy 
(STM) [4, 5]. They have been extensively studied [6], 
and play a key role in Kondo-Andreev screening [7] and 
mediating long-ranged antiferromagnetic exchange inter¬ 
actions [8]. Arrays of magnetic impurities with overlap¬ 
ping YSR states have recently attracted attention as a 
way to engineer a topological superconductor with Ma¬ 
jorana edge modes [9-16]. Realizing a Majorana mode is 
not only of fundamental interest, but also of technologi¬ 
cal importance, as their non-Abelian braiding statistics is 
key to promising proposals for a quantum computer [17]. 

Much excitement has therefore been generated by re¬ 
cent claims of the observation of Majorana modes in a 
ferromagnetic wire on a superconducting substrate [18]. 
Specifically, STM measurements reveal a sharply local¬ 
ized peak at zero bias in the local density of states 
(LDOS) at the ends of atomically thin iron wires grown 
on a lead surface, which are interpreted as Majorana 
bound states. Several theoretical calculations [19-21] 
have shown that Majorana inodes could occur in this 
setting, as had earlier been suggested in mesoscopic ge¬ 
ometries [22-26]. 

A characteristic aspect of the experimental results is 
that the LDOS in the wire shows pronounced subgap 
features [18], which resemble the spectrum of YSR states 
observed in isolated magnetic impurities or dimers [4, 5]. 
This motivates us to go beyond models of a spherical 
magnetic impurity [1-3, 6, 27], and consider the ferro¬ 
magnetic wire itself as a highly anisotropic scattering 
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Figure 1. (Color online) A wire-like magnetic impurity is em¬ 
bedded in a bulk superconductor with gap A. The wire has 
length L and radius a , with half-ellipsoidal caps at each end. 
We schematically show a YSR state as the scattering of an 
electron (solid line) at the impurity followed by its Andreev 
reflection as a hole (dashed line), and its subsequent scat¬ 
tering and Andreev reflection. Electrons and holes undergo 
conventional reflection at the magnetic impurity, described by 
the scattering matrices S and S, respectively. The difference 
in phase for electron- and hole-scattering off the impurity is 
compensated by the energy-dependent phase acquired by An¬ 
dreev reflection in the superconductor. 


center. By focusing upon the physics at subgap energy 
scales, we avoid the complexity of quasi-two-dimensional 
tight-binding models of the wire-like impurity [18, 21], 
yielding a transparent and intuitive understanding of the 
YSR states in this system. In particular, we wish to 
investigate whether or not the features interpreted as 
Majorana modes in Ref. 18 could have a conventional 
nontopological origin. 

In this letter we present a study of the YSR states 
in a ferromagnetic wire-like impurity with smooth ends, 
shown schematically in Fig. (1). We start by developing 
a general approach for obtaining the energies and wave- 
functions of the YSR states in terms of an eigenvalue 
problem involving the normal-state scattering matrices. 
We apply this formalism to a long wire-like ferromagnetic 
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impurity and calculate the energy spectra analytically us¬ 
ing a semiclassical approximation. This reveals that the 
YSR bands in the wire typically possess van Hove sin¬ 
gularities (VHSs) which appear as peaks in the LDOS 
of the wire. As the wire tapers to its end, these fea¬ 
tures move towards the gap edges. For sufficiently strong 
magnetization of the wire, the VHSs of the two particle- 
hole conjugated branches cross at zero energy, giving a 
sharply-localized peak in the LDOS near the wire end. 
We confirm this scenario by numerical computation of 
the scattering matrix of the impurity. The numerical re¬ 
sults for the LDOS qualitatively resemble key features 
of the experiment [18]. Finally, we show that including 
spin-orbit coupling in the impurity allows the appearance 
of Majorana modes inside a small energy gap. 

Scattering matrix formalism for bound states. We 
first examine a bound state wavefunction localized at a 
magnetic impurity in a superconductor, with energy ly¬ 
ing within the gap, i.e. \E\ < A [28]. Parametrizing the 
bound state energy through E = A cos <p, away from the 
impurity the wavefunction satisfies 


h 2 \/ 2 

2m 



if = A {—iT y + cos pT z )if , 


( 1 ) 


where m is the effective mass and /.i the chemical po¬ 
tential in the bulk superconductor. The solution can 
be expanded in eigenstates (l,e ±lv ) T of the matrix 
(—iTy + cos pr z ) and angular momentum channels 

V’(r) = X! { h< i~ ] (k-r)ai, m><T f ^ 

l,m,a ' ' 

+h\ + \k + r)b l ^ l7 ( ) | Yr(r )(a , (2) 

where a; jTnjtT and b bmi<7 are the scattering coefficients, ( a 
is a spinor, and the Hankel functions obey 

^ 2 ——f pj h\ ± \k±r) = ±iAsuuph\ ± \k±r ), (3) 



with k± = [2m(fi ± iAsiiup)} 1 / 2 /h. In the asymptotic 
limit \k±\r 1, and with /. i A, we approximate 


if(r) 



E 

Z,ra,cr 


i p -i[kFr-(l+ 1 ) tt / 2 ] 

^ ° Ui L,m,<7 



+el [k F r-(l+l),/2] bima 


1 


Yi m (i) (4) 


where vp = hkpfm is the Fermi velocity and kp = 
y/2m/a/h is the Fermi wavevector. 

For r <C £ = hvp/A we may ignore the exponential 
decay in Eq. (4), and the form of if(r) in the electron 
sector is then identical to the wavefunction for scat¬ 
tering off the impurity in the normal state. The co¬ 
efficients obey b = Sa, where we have introduced a 
matrix-vector notation for the normal-state scattering 
matrix S = and the vectors a = a brn ^ 


and b = b brrlja . Similarly, the hole channel wavefunction 
is a normal scattering state for the time-reversed Hamil¬ 
tonian H = a y H*a y . Therefore, the vectors a, b also 
satisfy e~ 2lip b = Sa , where the time-reversed scattering 
matrix S is given by S = Aa y S T a y A^ and A is the matrix 
that time-reverses the orbital indices. We hence obtain 
a and the phase p, which determines the bound state 
energy E = A cos p by solving the eigenvalue problem 

e~ 2itp a = S^Sa. (5) 

The other scattering amplitude b is obtained from a using 
b = Sa. Equation (5) generalizes the theory of YSR 
bound states [1-3] to arbitrary geometries. 

Semiclassical theory for a ferromagnetic wire. As an 
application of Eq. (5), we turn our attention to a wire¬ 
like impurity, which we set to lie along the 2 -axis. The 
scattering matrix can be estimated semiclassically, as¬ 
suming that the wire is smooth on the scale of kf, 1 . This 
is satisfied for the experiment reported in [18], where in 
the lead substrate we have kp 1 ss 0.13apb, with apb the 
lattice constant [29]. The small value of k F 1 makes it rea¬ 
sonable to assume that the scattering potential changes 
sufficiently slowly along the wire such that the paral¬ 
lel component of quasiparticle momentum k z is approx¬ 
imately conserved. For a cylindrically symmetric wire, 
the calculation of the bound-state energies for fixed k z 
reduces to the bound-state problem for an isotropic im¬ 
purity in two dimensions. Setting the magnetic moment 
along the 2 -axis, the scattering matrix in azimuthal an¬ 
gular momentum channel to is 

S m (k z , Z) = ; ( 6 ) 

where <p m (k z ,z) [ab m (k z ,z)] is the spin-independent 
(spin-dependent) component of the phase shift acquired 
by a scattered spin-er electron, and a is positive (nega¬ 
tive) according as the spin is aligned (antialigned) with 
the impurity magnetic moment. In the supplemental ma¬ 
terial we show that 4> m and b m satisfy [30] 

tan (4> m + ab m ) 

_ \k(j\Jm(k_[_a)J rn {k a of) k±^ </ m {k j_o.) (k a a) 

I k<j | b m (fc_L a)J m (k a a) k±A m (kxaf)J m {k a of) 

where k= kp — k 2 , k 2 = k\ — k 2 F V a l\i with V a the scat¬ 
tering potential experienced by spin-cr particles, J m (x) 
and Y m (x) are Bessel functions, and J m is J m (cc) or 
the modified Bessel function J m (|a;|) according as k 2 is 
positive or negative, respectively. The radius a of the 
wire vanishes towards the ends, giving the 2 -dependence 
in Eq. (6). Reflection and spin-rotation invariance imply 
4> m (-k z ,-z) = (fm(k z ,z) and b m (-k z ,-z ) = b m (k z ,z). 
The time-reversed scattering matrix is 

Sm(k z , 2) = e i<t>m{k z ,z)e-i b m (k„ z )<7 z _ (g) 

The eigenstates of S^Sm = exp(—2ib m a z ) yield the en¬ 
ergy of the spin-cr state in the to channel 

Em,a(k z ,z) = crAcos [b m (k z ,z)]. (9) 
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Figure 2. (Color online) Semiclassical results for the m = 0 
YSR bands in a magnetic wire, (a) Dispersion of the bands 
for different wire thicknesses a. (b) Dependence of the bound 
state energy on the wire thickness a at various fixed k z . The 
spin-']" (f) branch is represented by the thick (thin) lines. We 
set the potentials Vf = 0.04 fi and If = 0.6/r. 


We consider a magnetic potential sufficient to cause 
the m = 0 YSR bands to cross zero energy. Typi¬ 
cal results for the m = 0 bands are shown in Fig. (2); 
results for to > 0 are provided in the supplemental 
material [30]. The YSR bands generically have VHSs, 
which produce strong peaks in the DOS [16] as shown 
in Fig. (3)(a). In particular, the peak at ~ —0.5A cor¬ 
responds to the k z = 0 VHS, while the singularity at 
~ 0.8A matches the feature centered about k z /kF ~ 0.7. 
These VHSs are quite robust to changes in the scatter¬ 
ing potentials [30]. The reduction of the radius at the 
ends of the wire weakens the effective strength of the im¬ 
purity, causing the bands to shift towards the gap edge, 
as shown in Fig. (2)(b). The VHSs will therefore also 
shift, and so the LDOS near the wire ends is qualita¬ 
tively different from that in the middle. We estimate 
the length scale of this reconstruction by modeling the 
YSR band near the VHS with the effective Hamiltonian 
H = —h 2 d 2 /2fh + (z — zq)Vo, where the linear potential 
accounts for the movement of the VHS, which crosses 
zero energy at 2 = z 0 - The zero-energy solution is the 
Airy function Ai([z — zo]/4ff), with characteristic length 
scale l e g = {ti 2 /2inV 0 ) 1/3 . Estimating the effective mass 
to ~ kp/A and the rate of change in the position of the 
VHS Vo ~ Akp, we deduce that the peak in the LDOS 
will be localized over a region comparable to Z e ff ~ kp 1 . 
The peak at zero energy is expected to be particularly 
pronounced, due to the constructive effect of the cross¬ 
ing of the VHSs of the spin-f and -f branches. 

Numerical results: Here we confirm the intuitive 
semiclassical picture of the YSR physics by numerical 
computation of the scattering matrix, and furthermore 
calculate the LDOS along the wire. Our starting point 
is the relationship between the scattering matrix S and 
the T matrix [31] 

r = 5i t i> - 47rz k F (l, m\T a \l r , to) , (10) 

where (r|Z,m) = ji{kFr)Y l m (r)/y/2n. The T-matrix 
for spin-er scattering is formally given by T a = V a [l — 
GoV j]- 1 , where G 0 (r,r') = — e lkF ^ T ~ r '^ / ^ii\r — r'| is the 
free Green’s function. The T matrix is evaluated by 



Figure 3. (Color online) Numerical results for the m = 0 YSR 
states at the wire-like impurity, (a) Comparison of the total 
spin-f subgap DOS for the infinite wire (“bulk”) and the finite 
impurity, (b) Energy-resolved LDOS as a function of position 
z along the wire, with maximum value normalized to one. (c) 
Zero-energy LDOS as a function of position along the wire. 
We take the same potential for the wire as in Fig. (2). The 
wire length is chosen to be L = 20 a, and the Fermi wavevector 
is related to the radius by kpa = 6. We restrict the expansion 
of the T matrix to angular momentum channels with l < 70. 


expanding the Green’s function and potential in terms 
of angular momentum eigenstates on a real-space grid 
inside the wire. Having calculated the scattering ma¬ 
trix, we then use Eq. (5) to obtain the energy levels, 
from which we extract the total DOS of the wire plot¬ 
ted in Fig. (3) (a). There is excellent agreement with the 
thermally-broadened bulk semiclassical results. 

To obtain the local DOS we compute the wavefunctions 
within the impurity. In the T-matrix formalism, the scat¬ 
tering wavefunction \ip) is related to the bare wavefunc- 
tion \<j)) by | ip) = (1 — GoH) _1 |0). We obtain the bound 
states by choosing (p n - m ,a{r) = J2i a n -,i,m,a{r\l,m), where 
a n is the nth eigenvector of Eq. (5) corresponding to a 
state with energy E n . ma = aE n . m . The particle- and 
hole-components of the spin-f bound state wavefunctions 
are then 

Un-m, t(r) = y^(r|(l - , (Ha) 

l 

Vn-m, t(r) = ^(r|(l - G 0 Vi)~ 1 \l,m)a n .i^ . (lib) 
i 

Summing over the bound states and integrating over the 
direction perpendicular to the wire axis, we obtain the 
LDOS as a function of the position along the wire 

p m (z,E) oc / d 2 r ± ^2{\u„. m ^(r)\ 2 5(E - E n . m ) 

Jr ± <a(z) n 

+ k ;m , t (r)| 2 S(E + E n;m )} jE ^ - EIJ 12) 

The factor of yjE 2 — E%. m is due to the greater delocal¬ 
ization of the bound state near the gap edges. In eval- 
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uating the LDOS, we include thermal broadening of the 
^-functions at a temperature ksT = 0.1A. 

We plot the LDOS in Fig. (3)(b). In the bulk of the 
wire (i.e. 2 a < z < 18a) the LDOS shows very weak z- 
dependence. The main features are broad peaks centered 
at \E\ ~ 0.5A and 0.8A, which correspond to the VHS 
of the wire-limit bands. Within the end caps (i.e. z < 
2 a, z > 18a), however, there is substantial modification 
of the LDOS: The peaks at ±0.5A move towards =fA, 
and their crossing at zero energy gives a strong peak in 
the LDOS, as anticipated above. The zero-energy LDOS 
along the wire is shown in Fig. (3)(c), and reveals that 
this peak is sharply localized at the ends. The LDOS 
for the m =1,2 channels show similar behaviour. For 
to > 2, however, the potential is not strong enough to 
cause the bulk YSR bands to cross zero energy. As such, 
the VHSs of these bands do not shift through zero energy 
near the end of the wire, and they do not contribute to 
the enhancement of the zero-energy LDOS there [30]. 

Our model results show strong qualitative similarities 
to the experimental data in Ref. 18. Specifically, we re¬ 
produce the subgap peaks in the bulk of the wire and the 
sharply localized zero-bias peaks at the wire ends. In our 
model, however, the latter feature does not arise from a 
bound state (Majorana or otherwise) distinct from the 
bulk YSR bands, but rather originates from the shifting 
of the bulk bands at the ends of the wire. The result¬ 
ing zero-energy peak in the LDOS at the wire end is 
rather generic, requiring only that the impurity poten¬ 
tial be strong enough that the VHSs of the YSR bands 
cross zero energy as the radius of the wire vanishes. A 
quantitative comparison with [18] would require us to 
determine the relative contributions of each m to the 
STM measurements, which is well beyond the scope of 
the current work. Nevertheless, our considerations show 
that the observed zero-energy LDOS features at the wire 
ends cannot be regarded as an unambiguous signature of 
a Majorana mode. 

Majorana modes in spin-orbit coupled wires: So far 
we have restricted our attention to purely ferromagnetic 
wire-like impurities. Realizing a topologically non-trivial 
state in such a wire, however, requires SOC either within 
the wire [32] or the superconductor [22, 24]. Here we 
modify our semiclassical formalism to account for the 
former scenario. Neglecting the ^-dependence of the wire 
profile, we generalize the scattering matrix Eq. (6) and its 
time-reversed form Eq. (8) as S m (k z ) = exp (i[</> m (k z ) + 
b m (/c z )-er]) and S m (k z ) = exp(i[<f> m (-k z )-b m {-k z )-(r]), 
respectively, where the orientation of b m (k z ) depends 
on k z due to the SOC in the impurity. For a mirror- 
symmetric geometry, we choose (f>m{k z ) = 4> m (—k z ) and 
b m (k z ) = b m (k z ) sm9 m (k z )x + b m (k z ) cos9 m (k z )z where 


b m {k z ) = b m {—k z ) and 9 m (k z ) = - 9 m (-k z ) is the angle 
of the momentum-dependent magnetic moment with the 
wire axis. Solving Eq. (5), we obtain the YSR bands 

E m ,±{k z ) = ±A^ 1 - cos 1 2 9 m [k z ) sin 2 b m (k z ), (13) 

showing that the SOC gaps the spectrum Eq. (9) of the 
wire. 

The wire with SOC is in symmetry class D 1 which is 
characterized by a Z 2 topological number, specifically a 
Pfaffian. In our semiclassical model this is determined by 
the spectrum at k z = 0. Since the SOC is irrelevant here, 
the Pfaffian W = sgn{U 0 ,t(fe = 0)}, where E m ^ a (k z ) 
is the dispersion in the absence of SOC [Eq. (9)]. The 
topologically non-trivial state corresponds to W = — 1. 
From Fig. (2) we see that for the parameters chosen 
in our numerical study, adding SOC places the wire in 
the topological regime, allowing it to support Majorana 
modes at its ends. These zero-energy Majorana modes 
would appear in the background of the low-energy LDOS 
in Fig. (3)(b). At finite temperatures, however, all these 
subgap features will broaden considerably, making exper¬ 
imental interpretation difficult [20]. 

Conclusions. In this letter we have studied the YSR 
states in a wire-like ferromagnetic impurity embedded in 
a superconductor. Such systems are currently of great 
interest as a possible platform for realizing Majorana 
modes [18-21]. Using both analytical and numerical 
methods, we find that the YSR bands typically display 
VHSs, which appear as pronounced peaks in the wire 
LDOS, and which move towards the gap edges at the 
ends of the wire. For sufficiently strong magnetic scat¬ 
tering, the two YSR branches cross near the wire ends, 
giving rise to sharply-localized features in the LDOS at 
zero energy. Addition of SOC in the wire can stabilize 
a topological regime, but any Majorana features would 
appear in the background of an already elevated LDOS 
near the wire ends. The localized zero-energy peak in 
the LDOS obtained here is qualitatively consistent with 
recent STM measurements [18] and is independent of the 
existence of Majorana modes. A conclusive demonstra¬ 
tion of Majorana modes will need to distinguish them 
from the zero-energy conductance feature presented here. 
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I. SEMICLASSICAL THEORY OF THE INFINITE FERROMAGNETIC WIRE 

A. Scattering phase shifts 

We consider an infinite ferromagnetic wire impurity of radius a , lying along the z axis. The Hamiltonian for spin-u 
particles is written 


/i 2 V 2 

H=~— - p + V a Q(a-p) (SI) 

2m 

where p = K 2 k 2 F /2m is the chemical potential and V a is the impurity potential of the wire. Due to rotational symmetry 
about the z-axis and translational symmetry along the wire, solutions of Eq. (SI) have the form 

W(r )='l>m,k z ApV m *e ikzZ , (S2) 


where m is the azimuthal quantum number and k z is the momentum along the 2 axis. A general ansatz for the radial 
function ip m ,k z ,cr{p) is 

1^m,k z ,(7{p) ~ n) {Qm,k z ,aJm[kA-P) T fim,k z ,a J m (^1/^)} 

+0(a - p){Q(h 2 k\/2m - V a )c m ^ z ^J m {\k a \p) + <d(-H 2 k 2 ± /2m + V a )c m ,k z ,alm(\ka\p)} ■ (S3) 


Here J m (x) and Y m (x ) are Bessel functions of the first and second kinds, respectively, while I m (x) is the modified 
Bessel function of the first kind. The wavevector outside of the wire is k F = k 2 F — k 2 : while k 2 = k\ — 2mV a /h 2 is 
the wavevector or decay length within the cylinder, according as the term inside the absolute value sign is positive or 
negative, respectively. The coefficients a m ^k zF , Pm,k zi a, and c m ^k z ,a are C numbers. 

Far away from the scattering center (i.e. p^> a) we have the asymptotic form 


i>m,k z ,a{p) ~ a m ,k z ,a^== cos(k±p - TO7t/2 - 7t/4) + p m ,k z ,a^, sin(fc_i_p - 77171-/2 - 7 t/4) 


V k ±p 


V k ±p 


oc 


V k ±p 


cos(fcj_p - rmr/2 - it/A + S m>kzta ) 


(S4) 


where S m ,k z ,<j = <t>m(k z )+<xb m (k z ) = - arctan (P m ,k z ,a/ a m,k z ,a) is the scattering phase shift. The ratio (3 m ,k z ,a/am,k z ,a 
is determined from the boundary conditions at the edge of the wire 

i>m,k z ,a{p = a + ) = 1pm,k z ,<r(p = a ~ ) > (S5) 

dp^m,k z ,a{p = a + ) = d p l/j m ,k z ,a(P = a~) . (S6) 


After straightforward algebra, we deduce 


tan {S m ,k z ,a) = tan (</> m (fc 2 ) + ab m (k z )) 


This is given by Eq. (7) in the main text. 


ka'}n, (kci) J rl (kfj u) kJ n (kci)J n (k(jCi^ 
kfjYji {k(i) J' n {kfj(x) kY n {kct) J n (/c<ju) 

I kcr | Jn (kci)I n (| k a |u) kJ n {kci)I n {\k(j |u) 

\K\Y n (ka)r n {\k a \a) ~ kY^{ka)I n {\k a \a) ’ 


h 2 k\/2m > V a 
h 2 k\/2m < VQ 


(S7) 


B. YSR bands 

The YSR bands are obtained from the formula [Eq. (9) of main text] 


E m ,o{k z ,z) = dAcos[& m (fc z , 2 )] 


(S8) 
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Figure SI. Semiclassical results for the dispersion of 0 < m < 5 spin-L YSR bands in a magnetic wire for different wire 
thicknesses a. We set the potentials Vf = 0.04 p and Vj_ = 0.6 /n. The m = 0 results are given as Fig. (2)(a) of the main text. 



Figure S2. Semiclassical results for the dispersion of the m = 0 spin-')' YSR band upon varying (a) the nonmagnetic potential 
U and (b) the magnetic potential V. We set fea = 6 in all figures. In (a) we take V = 0.28/r while in (b) we have U = 0.32/r. 


where b m (k z ,z ) is calculated from Eq. (S7). We plot the 0 < m < 5 spin-f bands as a function of k z for different 
wire radius a in Fig. (SI). We note that all bands show a pronounced van Hove singularity (VHS) at k z = 0. The 
impurity potentials are such that the m < 2 bands close the gap in the bulk of the wire, i.e. they cross zero energy, 
with the VHS located at negative energies. As the width of the wire is reduced, these channels are therefore expected 
to contribute to the enhancement of the local density of states (LDOS) at zero energy near the end of the wire. On 
the other hand, the impurity is not strong enough to cause the m > 2 YSR bands to cross zero energy, and these 
channels therefore do not contribute to the enhancement of the zero-energy LDOS at the ends of the wire. 

The qualitative form of the YSR bands is robust to significant changes in the nonmagnetic component of the 
potential U = + Vj_) as shown in Fig. (S2)(a): increasing U by 50% produces only minor changes in the location 

of the van Hove singularities. The bands are more sensitive to changes in the magnetic potential V = — V|) as 

seen in Fig. (S2)(b), but still there is a large range of V values for which the important van Hove singularity at k z = 0 
is present. This indicates the generic nature of our results and the absence of fine-tuning. 
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Figure S3. Numerical results for the energy-resolved LDOS along the wire for the different m channels: (a) m = 0 [also Fig. 
(3)(b) of the main text]; (b) m = 1; (c) m = 2; (d) m = 3; (e) m = 4; and (f) m = 5. In each figure the LDOS is normalized 
by its largest value. The potentials in the ferromagnetic wire are Vf = 0.04/i and Vj. = 0.6 fj,. 



II. NUMERICAL RESULTS FOR WIRE LDOS 


The subgap LDOS in the wire for the to = 0... 5 YSR states is shown in Fig. (S3). As for the to = 0 results shown 
in the main text as Fig. (3)(b) [and reproduced here in Fig. (S3) (a)], we observe a significant difference between the 
states in the bulk and at the wire ends: the bulk is characterized by broad peaks centered at the location of the VHSs 
observed in Fig. (SI), but within the rounded ends of the wire these peaks move towards the gap edges. For the to < 2 
states, the k z = 0 VHS is located below zero energy, and so the shifting of the bands leads to an enhancement of the 
zero-energy LDOS at the wire ends. The potential is not strong to cause the to > 2 YSR bands to cross zero energy, 
however, and so the VHSs are not shifted through zero energy at the wire ends and these bands do not contribute 
to the zero-energy LDOS enhancement. Nevertheless, these bands would naturally explain the subgap features at 
nonzero energy observed in addition to the zero-energy peak at the wire ends in [18]. However, to quantitatively 
model these experiments we would need to be able to weight the contributions from the different to channels. This 
requires a detailed understanding of the coupling to the STM tip which is well beyond the scope of the current work. 













